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FIELD OF THE INVENTION 

The present invention relates to a method for estimating a propagation channel in the 
presence of transmit beamforming. It further relates to an estimator implementing said 
method. 

5 Said method is suited to be used in particular within the UMTS standard and 

especially within a conventional RAKE receiver. 

BACKGROUND OF THE INVENTION 

A UMTS communication system, defined in the UMTS standard (see the "3 GPP 

10 Technical Specifications, Rel. 99, http://www.3gpp.org/specs/specs.html" ^ incorporated 
herein by regerence, comprises at least a base station (BS) and a mobile phone terminal, 
referred to as user equipment (UE) throughout the document. Both the BS and the UE 
includes a transmitter and a receiver. The BS sends signals by means of the BS transmitter 
to the UE receiver in the downlink. The UE sends signals by means of the UE transmitter to 

15 the BS receiver in the uplink. In uplink and downlink communication both the UE and BS 
map the data to be transmitted into logical channels. 

In particular the UMTS standard provides two logical channels for the downlink 
communication, a first channel denoted common pilot channel (CPICH) consisting of a 
continuous sequence of known pilot symbols and a second logical channel denoted 

20 dedicated physical channel (DPCH). Said second logical channel DPCH consists of a first 
logical sub-channel denoted dedicated physical data channel (DPDCH), that carries the 
information data symbols intended for the user concerned, time multiplexed with a second 
logical sub-channel denoted dedicated physical control channel (DPCCH), that carries 
known user-dedicated pilot symbols. 

25 Transmit beamforming can be employed at the BS during downlink communication to 

improve the reception at the UE of the user information data mapped into the logical 
channel DPDCH. To perform transmit beamforming, the BS needs to be equipped with 
multiple antennas. In this way it is possible to provide optimal spatial filtering by shaping 
beams in the directions of the user concerned while setting nulls in other users directions to 
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mitigate the interference to other users. Hence the use of transmit beamforming can 
significantly improve the efficiency in the use of the spectrum, allowing a larger number of 
users to share the same spectrum, thanks to an improved user-dedicated data reception and 
reduced interference. While transmit beamforming is to be employed for the transmission of 
5 the logical channel DPCH, it is not suited for the logical channel CPICH transmission that is 
to be uniformly broadcast in time and space (i.e. in all directions) to all users. As seen from 
the UE receiver, the transmission of the logical channel DPCH occurs through a so called 
DPCH propagation channel consisting of the cascade of the BS transmitter chain, including 
the beamforming antenna system, the actual propagation channel from the BS to the UE and 

10 the UE receiver chain. Unlike the logical channel CPICH, downlink transmission occurs 
through a so-called CPICH propagation channel consisting of the cascade of the BS 
transmitter chain, including an omnidirectional antenna system in place of the beamforming 
antenna system, of the actual propagation channel from the BS to the UE and the receiver 
chain. For the sake of brevity we shall refer to a DPCH channel and a CPICH channel and 

15 the above DPCH propagation channel and CPICH propagation channel respectively 

throughout the present document. The UE receiver needs to estimate the DPCH channel in 
order to demodulate, decode, and recover the information data contained in the DPDCH. In 
the absence of transmit beamfonning the CPICH channel and the DPCH channel are the 
same (up to an uninfluential real scalar factor due to the different transmit power allocated 

20 to the DPDCH, DPCCH and CPICH logical channels). The continuous provision of known 
pilot symbols of the CPICH logical channel can be easily used at the UE receiver to 
estimate the CPICH channel, which is the same as the DPCH channel, and then use this 
estimate for DPDCH information data recovery. Unfortunately, when transmit beamforming 
is applied, the known pilots of the CPICH logical channel yet provide estimates for the 

25 CPICH channel, but those estimates cannot be directly applied as channel estimates for the 
DPCH channel that is different from the CPICH channel. Then the UE receiver has to rely 
on different techniques to estimate the DPDCH channel in order to recover the information 
data. 

30 SUMMARY OF THE INVENTION 

Accordingly, it is an object of the invention to provide a method for DPDCH channel 
estimation in the presence of transmit beamforming optimally exploiting all information 
available from the logical channels DPCH and CPICH, while limiting computational 
complexity, and an associated estimator implementing said method. 
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To this end, according to a first object of the invention, there is provided a method for 
estimating a propagation channel in the presence of transmit beamforming as claimed in 
claim 1. 

5 In addition, according to a second object of the invention, there is provided an 

estimator for estimating a propagation channel in the presence of transmit beamforming as 
claimed in claim 15. 

In a first non-limitative embodiment of the invention, the method is characterized as 
claimed in claim 2. 

1° In a second non-limitative embodiment of the invention, the method is characterized 

as claimed in claim 3. 

In a first non-limitative embodiment of the invention, the estimator is characterized 
as claimed in claim 16. 

In a second non-limitative embodiment of the invention, the estimator is 
15 characterized as claimed in claim 17. 

As we will see in detail further on, this method enables the UE receiver to cope with 
transmit beamforming applied to DPCH transmission by jointly exploiting the CPICH and 
DPCH logical channels information available at the UE. The first embodiment of the 

20 method optimally exploits this information according to the maximum-a-posteriori criterion, 
which is well known to the person skilled in the art, to estimate the DPDCH propagation 
channel multipath coefficients as described above. We shall show that in the presence of 
transmit beamforming, each multipath component coefficient of the DPCH channel differs 
from the corresponding multipath component coefficient of the CPICH channel by an 

25 unknown complex scale factor already referred to as beamforming complex factor. The 
second embodiment of the method provides an estimate of that beamforming complex 
factor. Then an estimate of the multipath component coefficient of the DPDCH propagation 
channel is built by multiplying each multipath component coefficient estimate of the CPICH 
channel by the related complex scale factor estimate. In this way a considerable computation 

30 saving is possible with respect to the first method. Moreover, the use of the CPICH channel 
estimates also allows to continuously track the channel variations with time due to the 
Doppler spread generated by the motion of the UE, thanks to the continuous provision of 
known pilot symbols of said CPICH logical channel. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention and additional features, which may be optionally used to implement 
the invention to advantage, are apparent from and elucidated with reference to the drawings 
described hereinafter. 

5 

Fig. 1 represents a diagram illustrating the main steps of said first method for DPDCH 
channel estimation according to the invention, 

Fig. 2 represents a diagram illustrating the main steps of said second method for 
DPDCH channel estimation according to the invention, 
10 - Fig. 3 illustrates an architecture schematic implementing the first method for DPDCH 
channel estimation according to the invention of Fig. 1, and 
Fig. 4 illustrates an architecture schematic implementing the second method for 
DPDCH channel estimation according to the invention of Fig. 2. 



15 DETAILED DESCRIPTION OF THE EMBODIMENTS 
1. Introduction 

In the following description, functions or constructions well-known to a person of 
ordinary skill in the art are not described in detail since they would unnecessarily obscure 
the invention. 

20 The present invention relates to a method and embodiments for the estimation of the 

propagation channel in a communication system when information data are transmitted 
through a beamforming system. Said methods are particularly suitable for implementation in 
a RAKE receiver within the UMTS standard when information data are transmitted from a 
base station BS to user equipment UE, here a mobile phone. 

25 The UMTS standard provides two logical channels for the downlink communication: a 

first channel denoted common pilot channel (CPICH) consisting of a continuous sequence 
of known pilot symbols and a second logical channel denoted dedicated physical channel 
(DPCH). Said second logical channel DPCH consists of a first logical sub-channel denoted 
dedicated physical data channel (DPDCH) that carries the information data intended for the 

30 user concerned, time- multiplexed with a second logical sub-channel denoted dedicated 

physical control channel (DPCCH) that carries known user-dedicated pilot symbols. Every 
time slot a DPCH is allocated, namely a DPDCH data field followed by DPCCH pilot field 
is allocated. 
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Transmit beamforming can be employed at the base station during downlink 
communication to improve the reception at the UE of the user information data mapped in 
the logical channel DPDCH. To perform transmit beamforming the BS needs to be equipped 
with multiple antennas. In this way it is possible to provide optimal spatial filtering by 
5 shaping beams in the directions of the user concerned while setting nulls in other users 

directions to mitigate the interference to other users. Hence the use of transmit beamforming 
can significantly improve the efficiency in the use of the spectrum, allowing a larger number 
of users to share the same spectrum, thanks to improved user-dedicated data reception and 
reduced interference. Whereas transmit beamforming is to be employed for the transmission 

10 of the logical channel DPCH, it is not suited for the logical channel CPICH transmission 

that is to be uniformly broadcast in time and space (i.e. in all directions) to all users. As seen 
from the UE receiver, the transmission of the logical channel DPCH occurs through a so- 
called DPCH channel consisting of the cascade of the BS transmitter chain, including the 
beamforming antenna system, the actual propagation channel from the BS to the UE and the 

15 UE receiver chain. In contrast, the logical channel CPICH is transmitted through a so-called • 
CPICH channel consisting of the cascade of the BS transmitter chain, including an 
omnidirectional antenna system in lieu of the beamforming antenna system, of the actual .t 
propagation channel from the BS to the UE and the receiver chain. The UE receiver needs to 
estimate the DPCH channel in order to demodulate, decode, and recover the information 

20 data contained in the DPDCH. Li the absence of transmit beamforming the CPICH channel 
and the DPCH channel are the same (up to a real scalar factor due to the different transmit 
power allocated to the DPDCH, DPCCH and CPICH logical channels). The continuous 
provision of known pilot symbols of the CPICH logical channel can easily be used at the 
UE receiver to estimate the CPICH channel, that is the same as the DPCH channel, and then 

25 use this estimate for DPDCH information data recovery. 

Unfortunately, when transmit beamforming is applied, the known pilots of the CPICH 
logical channel yet provide estimates for the CPICH channel, but these estimates cannot be 
directly applied as channel estimates for the DPCH channel because this is different from 
the CPICH channel. 

30 2. Fundamentals on Data Transmission and Reception in DS-CDMA 

In order for the information dedicated to a certain user to be transmitted? the source 
information data are properly mapped into a binary alphabet, i.e. into bits, via source- 
encoding techniques. Then these bits are channel-encoded by using suitable channel codes 
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in order to protect the information bits from all impairment causes (noise, interference) 
during the transmission. The channel-coded bits are then modulated by mapping bits into 
symbols according to a certain modulation scheme (e.g. QPSK modulation, well know by 
the person skilled in the art) where each symbol is represented by one or more bits (in 
5 QPSK modulation each symbol is represented by two bits). In DS-CDMA systems said 
symbols are then spread over a larger bandwidth by filtering with specific spreading 
sequences. In the literature spread symbols are referred to as chips. Notice that the duration 
of a chip period is shorter than the duration of a symbol period by a factor approximately 
equal to the bandwidth expansion factor of the spreading operation. Said factor is equal to 

10 the number of chips per symbol period and is referred to as spreading factor. Symbol 
sequences carrying information bits from different sources are spread with different 
spreading sequences in order for several users to share the same bandwidth and for the 
receivers to be able to discriminate and recover the desired user's symbol sequences based 
on the knowledge of the user specific spreading sequence. The obtained chips are then* . 

1 5 filtered by a pulse-shape filter (a root-raised cosine filter in the case of UMTS) and digital- 
to-analog (D/A) converted. The resulting analog signal is then modulated to radio frequency • 
and transmitted by the BS antenna(s). Similarly, when transmit beamforming is employed v 
the same sequence of chips filtered by the pulse-shape filter is sent to all the antennas after 
multiplication by an antenna-specific complex weight factor. The antenna-specific weight 

20 factors can be collected in a so-called beamforming weight vector. The choice of said 
beamforming weight vector determines the shape of the irradiation diagram, in terms of 
magnitude and phase, of the multiple-antenna system of the BS. The signals transmitted 
from the antennas undergo several scattering, diffraction, and reflection phenomena through 
the actual propagation channel from the BS antennas to the UB receiver antenna. Noise and 

25 interference from other users and neighboring BS further contribute to impairing the useful 
signal(s) for the UE concerned. The radio frequency signal received at the UE antenna is 
demodulated down to baseband (or intermediate frequency IF) and analog-to-digital (A/D) 
converted to generate a digital baseband signal. Then the UE receiver processes this 
baseband signal in order to recover the useful information data intended for the user 

30 concerned. For this purpose, the receiver needs to estimate the cascade of the transmitter 
chain, the actual wireless propagation channel and receiver chain up to the A/D conversion, 
associated with the transmission of the useful data. In the literature this cascade is in general 
simply referred to as the channel. In that which follows, a detailed mathematical description 
of the above operations is provided. 
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3. Signal and Channel Model 

General DS-CDMA signal and channel models are provided. Such models will then be 
applied to the specific cases of UMTS CPICH and DPCH signals and channels. 

3.1. General Signal Model 

5 We consider a general case of transmission of a signal s(t) through a multipath channel wit h 
impulse response h(t, r) . The continuous-time complex baseband received signal prior to 
A/D conversion at the receiver is modeled as: 

X0 = *(0+v(0 (1) 
where x(i) represents the part of the received signal that comprises the useful data and v(t) 
10 denotes the noise plus interference term. The signal x(t) is given by: 

*(*)= J/KM-rMr)dr (2) 
where h{t 9 r) represents a time- varying channel impulse response and 

z(0 = £a(*Mf-*r c ) (3) 

represents the transmitted signal after D/A conversion at the transmitter, withj c denoting 
1 5 the chip period and a(k) = s(\k/M\) d(k) , where M denotes the spreading factor, and [.J . . i - 1 
denotes the flooring operator. Moreover, *(») represents the w-th modulation symbol (e.g., : 
QPSK), d(k) denotes the Ar-th chip of the spreading sequence, and y^f) represents the pulse- 
shaping filter that limits the system bandwidth. Note that this model can be used by systems 
with several spreading layers, like DS-CDMA system defined by the standard UMTS, IMT- 
20 2000 and IS-95 (see "IMT-2000 http://ww.imt-2000.org/portal/index.asp^ incorporated 
herein by reference. 

The channel h(t,r) is assumed to follow the wide-sense stationary uncorrelated 
scattering model with Rayleigh fading and multipath response as described in "J. G. Proakis 
Digital Communications, NY: McGraw-Hill, 3 rd edition, 1995". 

25 /^r) = 2> p (0*(r-T,) (4) 

where P denotes the number of multipath components, and c p (t) and r p denotes the time- 
varying complex coefficient and the propagation delay associated with /?-th multipath 
component. The path complex coefficients c p (t) vary depending on the velocity of the UE 
and of the surrounding scattering objects (e.g., buildings, hills but also other moving 
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objects) with respect to the BS. When the BS uses Q antennas for transmit beamforming, the 
channel model expressed by (4) generalizes to 

W, r) = 2 w» a(0, )c p (t)S(r - r p ) (5) 

where 

5 - w = [w 0 , . . . , M/ e _! f denotes the transmit beamforming weight vector 

a(0) = [a 0 (0), . . . , denotes the antenna response vector in a direction 0 . 

The superscripts ( . ) T and ( . ) H denote transpose and Hermitian transpose respectively. 
Since in a practical implementation the vector w typically may vary at most at the slot rate 
for the sake of notation, its variations with time are not explicitly taken into account in 
10 equation (5). 

By inspecting equations (4) and (5) one may notice that, during a time slot, transmit 
beamforming only affects the amplitude of the path time-varying complex coefficients c p (t) 
by a complex factor w"a(0 p ) and not the delays r p . In other words, the channel's temporal 

structure is not affected by the antenna response vector, neither by the transmit 
15 beamforming. Thus the path time-varying complex coefficients j3 p Cp(t) = w H a(^ p )c p (t) and 

Cp(0 are the same in both case (4) and (5) up to a complex factor j8p = w H a(0 p ) 
representing the previously mentioned beamforming complex factor. Therefore, a channel 
estimate according to (5) can be built from estimates of the delays r p p of the coefficients 

Cp(t) and of P P for p=l,..., P. 

20 3.2. Received Signal Model for DPCH and CPICH 

We now apply the previous general channel model to the DPCH and CPICH channels 
considering the case of UMTS where the CPICH and DPCH channels are superposed. Then 
the continuous-time complex baseband received signal (1) can be written as 

^(0 = XcptcH (0 + * dpc * (0 + v(0 (6) 
25 where Xc P i C h (t) and Xd pC h (0 are the signal components associated with the CPICH and 
DPCH channels respectively, given by 
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where 

^(^) = fc p (0^-r p ) 

and 
where 

" (*) = WL*/256j)^(*) , 

- Sapidity), Sdpdiif) denote the n-th modulation symbol associated with the CPICH and the 
z-th modulation symbol associated with the DPCH respectively, and d cpich {k\ d dpch (Jc) 
donote the fc-th chips of the CPICH and DPCH spreading sequences respectively. 

3.3. Remarks 

It will be noted that the spreading factor M may not be the same for the CPICH and DPCH 
logical channels. In the UMTS standard the CPICH has a fixed spreading factor equal to 
256, whereas the DPCH has a spreading factor that can take on values in the set {4, 16, 32, 
64, 128, 256,512}. 

Furthermore, the CPICH symbols are all the same, i.e., Scpi C h(«) = A(l+j) for any n, 
where A is a real amplitude factor taking into account the power that the BS assigns to the 
CPICH logical channel. The a-priori knowledge at the UE receiver of the symbols 
comprising the CPICH can be exploited to estimate the channel hcp ich as expressed by (4). In 
the absence of transmit beamforming on the DPCH hcpi ch = h^h, so that such an estimate 
based on the CPICH symbols would suffice for the receiver to recover the desired user's 
information symbols on the DPCH. On the contrary, if transmit beamforming is applied to 
the DPCH, the a-priori knowledge of the CPICH symbols yet allows to easily estimate the 
channel hcp ich as expressed by (4), but not the channel hdpch associated with the DPCH 
expressed by (5) which differs from the hcp ich due to the presence of the unknown 
coefficients p p . Hence other approaches need to be addressed to provide an estimate of the 
channel hd pC h necessary to recover the information symbols transmitted over the DPDCH. 
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We have implicitly assumed above a digital baseband receiver where baseband signals 
are processed after A/D conversion. In other words the signal expressed by (6) is A/D 
converted, prior to all baseband processing. In particular we consider the UE receiver to be a 
discrete-time RAKE receiver (see J. G. Proakis Digital Communications, NY: McGraw Hill, 
5 3 rd edition, 1995 incorporated herein by reference). We shall notice that the methods hereby 
addressed can be applied to all types of baseband receivers relying on channel estimation 
though. 

4. Channel Estimation 

Since the delays r p are not affected by the presence of transmit beamforming, we shall 

10 assume that they are already recovered by usual delay estimation/tracking techniques (e.g. 
delay lock loop, see J. G. Proakis Digital Communications, NY: McGraw-Hill, 
3 edition, 1995 incorporated herein by reference and exploiting the knowledge of the 
CPICH symbols) and we shall focus on the estimation of the channel path coefficients c p (t) 
and Pp. In order for those coefficients to be estimated we consider the signals at each finger 

15 of the RAKE receiver after despreading of both the CPICH and the DPCH. 

For the sake of simplicity we shall limit the treatment to the single path/finger case, 
the multi-finger case being a straightforward generalization of it. Thus, we will drop the 
subscript p to denote the time-varying complex coefficient Cp(t) liereinafter also referred 
to as path fading coefficient, or simply path coefficient. Moreover, for the sake of brevity in 

20 the following we shall refer to "channel estimation" and "channel estimate" as the channel 
path coefficient estimation and estimate respectively. 

Furthermore, in a first instance, we also limit the analysis to the case where the 
DPCH spreading factor Mis equal to that of the CPICH. We shall generalize the analysis to 
other spreading factors successively. 

25 In addition, we need to model the previously mentioned interference and noise term. 

In particular the multi-access interference arises from the interference among different users 
and the interferences due to the CPICH over the DPCH and/or due to different data flows 
dedicated to the same user, and the inter-path interference is due to the interference between 
the impulse responses associated with paths at different delays r. 

30 Assuming that the whole interference and noise to be modeled as a zero-mean 

complex circularly symmetric additive white Gaussian noise (AWGN), the CPICH and 
DPCH de-spread signals at a RAKE finger are respectively given by: 
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y^(n) = yf H m<n)s dpdl {n) + v dpeh (ri) ( } 

where v dpch (n) and Vcpi C h(n) denote the noise-plus interference terms associated with the 
DPCH and CPICH respectively, assumed to be independent, c(n) is the discrete-time 
representation at symbol rate of the path time-varying complex coefficient previously 
5 denoted c(t) or Cp(t) , s dpc h(n) and Scpi C h(n) denote the symbols of the DPCH and CPICH 
logical channels respectively, at a symbol time n. Under the assumption that the CPICH and 
DPCH channels may vary every symbol time «, and that the interference and the noise can 
be modeled as a zero mean AWGN, we then build a maximum likelihood (ML) static 
instantaneous channel estimate at every symbol time n as follows: 
1Q ^(n) = s% tch {n) ycpich {n) = c(w) + v^/i) 

c dpch {n) = s^Wy^Qi) = w"a(0)c(*) + v dpch (n) 

where, s*(n) denotes the complex conjugated of s(n), Vcpi ch (n) and v dpch (n) represent the 
AWGN associated with the CPICH and DPCH channels respectively, and where, without 
loss of generality, we assumed the symbols of both channels to be normalized such that 
lsd P ch(n)| = |scpi ch (n)| = 1. In order for the estimates (8) to be computed one needs to know 

15 the transmitted symbols S* dpC h(n) and S*cpi Ch (n) . The symbols S*cpi C h(n) are all known at 
the receiver. Within the DPCH only the dedicated pilot symbols S* dpC ch(n) comprising the 
DPCCH are known, but not the data symbols comprising the DPDCH. Therefore, if 
DPDCH symbols S* dpdc h(n) are to be used in (8) they first need to be estimated by, e.g., a 
decision-direct technique. The ML estimates (8) are therefore formed from the observations 

20 ycpjch and y dpC h at the output of a finger in a data aided (DA) fashion when there are a 

priori known pilot symbols, as for the cases of CPICH and DPCCH, and in a decision direct 
(DD) fashion in the presence of unknown data symbols as for the DPDCH. The DD 
approach requires the use of hard decisions taken on the DPDCH at the RAKE receiver 
output as pilot symbols. Therefore, from (8) we have 

25 c dpech (n) = s dpccb (n)y dpcch (n) = w H z(0)c(n) + v dpcch (n) (9) 

for the DPCCH case (DA) and 

= K&thWydpdchW = w"a(0)c(*) + (n) (10) 
for the DPDCH case (DD) where (») represents the estimate of the n-th DPDCH symbol 
built as a hard decision taken at the RAKE receiver output. Equations (8)-(10) represent the 

30 basis for the derivation of the DPDCH channel estimation algorithms addressed in the 
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following. In order to support the next derivation, further details are provided concerning 
the physical behavior and related models for the CPICH and DPCH channels. 

4.1. Physical Interpretation of CPICH and DPCH Fading Channel Models 

Due to the large bandwidth of the UMTS W-CDMA system the channel from the BS to a 
5 UE can be modeled as comprising a set of discrete multipath components. In particular it is 
possible to approximate the channel impulse response as a linear superposition of a finite 
number of independent multipath components as expressed by equations (4) and (5). In 
what precedes we have also referred to a multipath component simply as a channel path. A 
path is characterized by a delay and a complex coefficient When the UE (and/or any other 

10 scattering object) is in motion with respect to the BS, each path complex coefficient may 
vary with time. Thus the overall channel impulse response varies with time as well, as 
expressed by equation (4) and following. This phenomenon is typically referred to as 
Doppler spread, from the Doppler effect arising when in a communication either the 
transmitter or the receiver moves one with respect to another. The channel Doppler spread 

15 represents the rate at which the channel impulse response varies, or focusing on a channel 
path, the rate at which the path complex coefficient varies with time. The Doppler spread 
increases with the velocity of the UE with respect to the BS. From a physical standpoint to 
each path of the CPICH channel corresponds a path on the DPCH channel with the same 
delay r and the same Doppler spread but with different complex-path amplitudes, namely 

20 c(n) and p c(n) respectively, according to equations (4) and (5). Then as long as the 

factorp is constant, both CPICH and DPCH channel paths vary at the rate of c(n) , namely 
they share the same Doppler spread. Hence the two complex paths' amplitudes temporal 
variations are the same while they always differ only by the beamforming complex factor p 
that is determined by the scalar product of the beamforming weight vector w and the 

25 antenna array response a(0) corresponding to the direction 9 of that path. 

4.2. Optimal Joint MAP CPICH-DPCH Based DPDCH Channel Estimation 

The optimal joint CPICH-DPCH MAP DPDCH channel estimation problem is formulated 
as a direct extension of the simpler CPICH MAP channel estimation problem addressed in 
"J. Baltersee, G. Fock, P. Schultz-Rittich, and H. Meyr, 'Terformance analysis of phasor 
30 estimation algorithms for FDD-UMTS RAKE receiver," IEEE & h Symp. on Spread 

Spectrum Technologies and Applications, NJTT, NJ, September 2000" incorporated herein 
by reference. 
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The general formulation of the MAP optimization criterion corresponds to the following 
optimization problem 

c MA p (*) = argirux p ( c (k) | y , s) 

namely c mp (k) is the argument that maximizes the a posteriori probability density function 
5 p(c(k) | y , s) of c(k) given the vector y of the observations y^ (n) and («) expressed by (7) 
and the vector s of the pilot symbols s cplch (n) and s dpch (n) . 

We shall see that even though optimal MAP DPDCH channel estimation is generally too 
computationally demanding for a feasible implementation within a realistic UE, its 
derivation provides the proper guidelines for other reduced complexity estimator for the 

10 DPDCH channel yet according to certain optimality criteria. 

Under the assumption that the channel paths are a Gaussian complex random process, 
the maximum a posteriori probability (MAP) DPDCH channel path coefficient estimate at 
time instant k is given by the conditional by expected value c MAP (k) of (k) = c dpdch (&) given 
the observation column vector c ML comprising all the available ML instantaneous channel 

15 path coefficient estimates c dpch (m) withm s*k, and c^n) defined according to equations 
(8)-(10),i.e.: 

£ mp (k) = E{c ML (k)\c m } (11) 

Equation (11) can be simply explained by observing that the ML estimates c ML {k) according 
to equations (8)-(10) are still Gaussian complex random variables with mean given by the 

20 true DPDCH channel c dpdeh (k) = J3c(k) and variance given by the variance of the noise term 
v dpch (k) . In order to achieve a MAP estimate the a-posteriori probability density function of 
Cml (*) (i-e., the probability density function of c m (k) conditioned to c ML ) has to be 
maximized. It is well known that the maximum of a Gaussian probability density function is 
reached in corresponding to the mean value of the distribution. Henceforth taking the 

25 expectation equation (1 1) maximizes the conditional probability density function yielding a 
MAP channel estimate. It can be shown (see e.g. H. Meyr, M. Moeneclay, and S. Fetchel, 
"Digital Communication Receivers: Synchronization, channel estimation and signal 
processing", John Wiley & Sons, New York, 1998" incorporated herein by reference) that 
a MAP optimal channel estimate (1 1) can be built by low-pass filtering a set of 

30 instantaneous ML channel estimates (8) as follows: 

CMAp(k) = f H c UL (12) 
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where f = [f(N - 1) ... /(0)] r in (1 2) is the vector of the low pass filter coefficients, with 
N = K cpldl + K dpcb , K dpch and AT^ being the numbers of the available ML instantaneous 
channel estimates of the DPCH, c dpeh (m) , and of the CPICH, c cpich {n) 9 respectively. Under the 
previous assumption that the channel paths are Rayleigh fading the optimal filter 
5 f = [f(N - 1) ... /(0)] r yielding a MAP channel estimate is given by: 

f = R« , V (13) 
where R~ = Efc^cj^} represents the covariance matrix of the ML channel estimates, 
r^. = B{Q^c dpdch {k)} represents the vector of the cross-correlation between the ML channel 

estimates and the true DPDCH channel at time instant k. 
10 The optimal filter f given by equation (13) corresponds to the well-known Wiener filter, 

also known as a linear prediction filter. Depending on whether the vector c m includes only 
past, or only future, or both past and future ML instantaneous channel estimates of the form 
(8)-(10) with respect to the DPDCH channel c^(£) , equation (13) gives rise to a forward, 

a backward, or a forward-backward linear prediction filter respectively. 

1 5 The optimal linear prediction filter f can be easily computed also in the more general ; 

case of Rice fading channels as shown in great detail in H. Meyr, M. Moeneclay, and 
S. Fetchel, "Digital Communication Receivers: Synchronization, channel estimation and 
signal processing", John Wiley & Sons, New York, 1998" incorporated herein by reference. 
Despite the fact that theory clearly indicates a closed form expression for the optimal 

20 solution to the problem of MAP channel estimation, as demonstrated above, a practical 
implementation of such an optimal estimator poses several complexity issues which 
generally prevent its feasibility within a practical UE architecture. In the next section we 
shall address the implementation of a MAP channel estimator. To this end we shall rely on 
the concepts of Wiener filtering, linear prediction, as well as adaptive and LMS-Kalman 

25 filtering oriented to MAP channel estimation along with fading channel models, which are 
considered to be well-known by the person of ordinary skill in the art. 

4.2. L Implementation of the Optimal DPDCH MAP Channel Estimation 
In order to implement the optimal joint CPICH-DPCH MAP DPDCH channel estimator in a 
practical UE the optimum Wiener FIR filter coefficients f = [f(N - 1) ... /(0)] r are to be 
30 computed. For this computation to be possible the statistics of the observation c ML and of the 
true channel c dpdch (k) are needed at the UE. These quantities can be estimated assuming 
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structured channel and signal models, such as the Jakes' fading model (see: J. G. Proakis 
Digital Communications, NY: McGraw-Hill, 3 rd edition, 1995 incorporated herein by 
reference), and AWGN model for the noise and interference terms. Then estimates of the 
model parameters are needed, i.e., the Doppler spread (or the UE velocity relative to the 
5 BS), the signal power and the interference-plus noise power spectral density. To save 

computations, the UE may select the optimal Wiener FIR filter from a set of pre-computed 
filter coefficients based upon the estimates of the parameters of the statistic channel model. 
Unfortunately, even neglecting the complexity involved by the parameter estimation, this 
approach is unfeasible. Indeed the number of filter coefficients to be stored would be 

10 overwhelming compared to the memory storage capabilities of a real UE. A sub-optimal 
solution would consist of using a reduced set of Wiener filters that can cope with the 
maximum envisaged channel Doppler spread (normalized with respect to the symbol 
duration). Different filters would be selected based only upon the signal-to-noise-plus 
interference ratio (SINR) estimate and no longer on the Doppler spread estimate. Indeed, 

1 5 estimating the Doppler spectrum is generally much harder than estimating the SINR. 
Moreover, SINR estimates are needed anyhow to perform closed loop power control. In 
order to guarantee that no information loss occurs at the maximum Doppler spread for a 
given spreading factor, these filters need to be designed to cope with the widest envisaged 
Doppler bandwidth. This solution yields significantly suboptimal performance at low 

20 Doppler spread (i.e. low UE velocities) and/or low spreading factors, because the filter 

bandwidth is much larger than the one of the channel path process in these cases, leading to 
poor noise suppression. 

An alternative approach consists of building unstructured estimates of the observation 
covariance matrix R tt and the channel autocorrelation sequence r^. , without assuming a 

25 specific statistic model. Once such estimates are available the linear prediction filter f can 
be built according to equation (13) by replacing R~ and with their estimates. In practice 
it would be more convenient to directly estimate the inverse of R~ to avoid the complication 
of inverting the estimate of the R K . For this purposes several techniques are available in the 
literature. 

30 A more practical solution when there are of known pilots, i.e. over the DPCCH and/or 

CPICH, or reliable hard decisions taken via DD mechanism over the DPDCH are available, 
consists of estimating directly the Wiener filter coefficients in adaptive fashion, by 
sufficient training. Adaptive filtering inherently provides tracking of the short-term and 
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long-term channel variations. Tracking and adaptation of the prediction filter can be 
performed in a continuous fashion or on a block basis, e.g. by adapting the FIR prediction 
filter every slot. 

A lower complexity approach than FIR Wiener FIR linear filtering consists of using a 
5 first order IIR LMS-Kalman as linear predictor. The HR loop gain factor needs to be smaller 
than one for the HR filter to be stable. Under the assumption of a parametric channel 
statistic model it can be shown that the optimal loop gain factor yet depends on the noise- 
plus interference term and the channel path statistics. Moreover, the optimal loop gain 
cannot be rendered smaller than one for high Doppler spread and large spreading factors, as 

10 it is demonstrated in: J. Baltersee, G. Fock, P. Schultz-Rittich, and H. Meyr, 'Terformance 
analysis of phasor estimation algorithms for FDD-UMTS RAKE receiver", IEEE & h Symp. 
on Spread Spectrum Techonologies and Applications, NUT, NJ, September 2000 
incorporated herein by reference. Therefore, in order to be stable, the LMS-Kalman linear 
predictor can no longer be optimal in those cases. Nevertheless, the LMS-Kalman linear 

15 predictor is still an optimal solution for low UE mobility (Doppler spread smaller than 222 
Hz) and small spreading factor (less than 128). 

The above considerations allow estimating the amount of complexity involved by the 
computation of the optimal prediction filter. This complexity represents a major hurdle for a 
feasible implementation of the MAP channel estimator within a realistic UE, so that other 

20 sub-optimal but less computationally intensive solutions need to be addressed. The presence 
of the expectation operator in the MAP estimator expressed by equation (11) indicates that 
any sub-optimal estimator should low-pass filter the instantaneous ML channel estimates 
expressed by equations (8)-(10). The design and implementation of a sub-optimal low-pass 
filter for the ML channel estimates (8)-(10) in lieu of the optimal Wiener filter has to 

25 achieve a trade-off between complexity saving and performance degradation with respect to 
the optimal MAP estimator. Very simple low-pass filter structure can be shown to achieve 
nearly optimal performance in most practical scenarios with significantly reduced 
complexity. These considerations are all about the practical design of a channel estimation 
mechanism suitable for a feasible UE. In the light of the previous derivation we shall 

30 introduce a low complexity but still nearly optimal DPDCH channel estimation algorithm. 



4.3. Joint CPICH-DPCH Based Low Complexity DPDCH Channel Estimation 

In section 3.1 we have seen that the channel paths associated with the CPICH and DPCH 
share a common structure. In particular, it has been shown that both channels associated 
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with the CPICH and DPCH behave in the same way up to said beamforming complex factor 
J3 due to the transmit beamforming. The knowledge of said complex factor J3 for all channel 
paths allows building an estimate of the DPDCH channel from an estimate of the CPICH 
channel. This way the DPDCH estimation problem splits into two simpler estimation 
5 problems, one of estimating p and another of estimating the CPICH channel c(k) . Optimal 
MAP channel estimation can be considered to estimate the CPICH channel c{k) with 
significant complexity reduction with respect to the complete joint CPICH-DPCH MAP 
DPDCH channel estimation problem. The optimal MAP CPICH channel estimate at time 
instant k is given by the expectation of the ML instantaneous CPICH estimate c^Ofc) given 

10 all available ML instantaneous CPICH estimates of the form (») for n * k . The solution to 
this problem has the same form as (12) and (13), where the observation column vector i m is 
to be redefined so as to comprise just the CPICH and not the DPCH ML instantaneous 
channel estimates (8) — (10), and the cross-correlation vector in equation (13) is to be 
redefined as = E{c m c*(*)} . A detailed derivation of a MAP CPICH channel estimator is 

15 provided in said reference J. Baltersee, G. Fock, P. Schultz-Rittich, and H. Meyr, 

'Terfonnance analysis of phasor estimation algorithms for FDD-UMTS RAKE receiver", 
IEEE 6 th Symp. on Spread Spectrum Techonologies and Applications, NJIT, NJ, 
September 2000. Similar considerations as those mentioned in section 4.3 apply also in this 
case. If the UE still cannot afford the complexity of a MAP CPICH channel estimate lower 

20 complexity approach can be pursued by adopting sub-optimal low pass filtering in lieu of 
the optimal linear prediction filter. In the following we shall refer to (k) as an estimate of 
the CPICH channel c(k) regardless the low pass filter employed for the estimation. Along 
with an estimate (k) of the CPICH channel c(k) , the UE has to compute an estimate p of 
said complex factor /?, to build an estimate of the DPDCH channel as c dpdeh (k) = pc cptcIl (k) . To 

25 this end, we shall derive an optimal ML estimator for/? that can exploit both the DPCCH 
pilot symbols in DA fashion and/or the DPDCH data in a DD fashion, jointly with the pilot 
symbols provided by the CPICH according to equations (8)-(10). 
4.3 J. Maximum Likelihood Estimation of the Beamforming Complex Factor p 
In order to build an ML estimate for the factor p we need to build ML estimates of the 

30 cross-correlation between DPCH and CPICH instantaneous ML channel estimates, and of 
the autocorrelation of the CPICH instantaneous ML channel estimates. For this purpose we 
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shall rely on the expressions of the ML instantaneous DPCH and CPICH channel estimates 
(8) — (10). The theoretical expressions for the DPCH-CPICH ML channel estimates cross- 
correlation and CPICH ML channel estimate autocorrelation are: 

^(0 = E{c^(«)c^(n-/)} (14) 

5 and 

^(O^E^OOc^ (»-/)} (15) 
respectively, where a correlation lag / * Ois to be selected for noise suppression, as it will be 
shown in the following. From equation (8), we have: 
^('0 = c(w)+v^(») 
c dpch (") = w*a(0)c(ii)+ (n) 
10 where we denote fi = w H a(0) . Under the assumption that the noise terms (n) and 

v dpch (n)are independent zero-mean additive white Gaussian noises, expressions (14) and (15) 

for/*0yieldto: 

^(0 = E{c^(n)c^(n-/)} = 

= 0E{c{n)c\n -/)} + >aE{c(«)v^(/i -/)} + E{v^(»)c>-/)} + E{v^(n)v^(«./)} = 
= /£{c(n)c>-/)} 

and 

^(0 = E{c^(«)c^(«-/)} = 
15 = E{c(»)c> -/)} + E{c( W )v^ e/ , (*-/)} + E{ V(3rfcft (n)c> -/)} + B{v cpieh (n)v^ (*-/)} = 

= E{c(n)c*(*-/)} 
respectively. Thus the complex factor y? is readily given by: 

P=^T7K for any /^O (16) 

A practical implementation requires the expectation operator E{.} to be replaced by a 
feasible estimator yielding estimates $ dc (l) and ^(Z) for the cross- and autocorrelation (14) 
20 and (15) respectively, to build an estimate fi for p according to equation (16). A simple and 
effective estimate for </> dc (l) and ^(i) is given by the sample temporal average, i.e.: 

4(0=^— Z^w^o*-/) (i7) 

and 
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«U0 = - Y, £ c Pid ,W)c' cptch {n-T) (18) 

cpich n=m 

where n = m denotes the time instant which the averages (17) and (18) start from. Notice that 
for / > oit is needed z +1 . The estimates (17) and (18) are consistent estimates and 
commonly known as sample estimators. They are also known to be ML estimates of the 
cross and auto-correlation (14) and (15) under the assumption of Gaussian observations, as 
in our case (see the hereby incorporate reference "B. Porat, Digital processing of random 
signals, Englewood Cliffs, NJ, Prentice-Hall, 1994"). They are also known to be biased 
estimates for finite K dpch , and asymptotically unbiased, i.e. unbiased as K dpch goes to infinity. 
Moreover, it is well known that the ratio of two ML estimates is still a ML estimate. Thus 
an ML estimate for p is given by: 

^ = Mr (19) 

-p 2-, C cptch OOCq** (» - 0 

"-cptch n=m 

where the dependence on the lag / is due to the finite number of symbols K dp£h and 
available in practice. Notice that in the case where K epbA = the estimator (19) takes the 
form: 

1 m £^(n)c^(H-0 

^ (/)= T7K = ^ < 2 °) 



TCC ^ 



The estimator (20) has the advantage of avoiding the need for computing the ratio 
Kcpici, 1 K dpch 311(1 therefore has a reduced complexity with respect to the general estimator 
(19). In order to further limit the complexity of estimator (20) one might avoid the presence 
of a DD mechanism to build the ML instantaneous channel estimates over the DPDCH 
according to equation (10). In this case only the known pilot symbols comprising the 
DPCCH would be used to form DPCCH ML instantaneous channel estimates expressed by 
equation (9) and it would result in K dpch = K dpcch . 

Since equation (16) holds for any / * 0 , p m (/) expressed by (1 9) and (20) represents 
an ML estimate of p for any / * 0 . It is also well known that a linear combination of ML 
estimates is still an ML estimate. Thus in a preferred non-limitative embodiment, in order to 
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improve the estimate of p , one may consider the weighted sum of ML estimates of the 
cross and auto correlation at different correlation lag / , i.e.: 

*~ m l : s t C-^ (2D 



1=1 M=m 



where the coefficients a, and 6, are to be chosen according to the accuracy (expressed in 
terms of mean square error) of each ML estimates ^(0 and ^(Q (see: B. Porat, "Digital 
processing of random signals", Englewood Cliffs, NJ, Prentice-Hall, 1994 incorporated 
herein by reference). 

Moreover, it can be shown that the estimators (17) and (18) are asymptotically 
efficient only for the first correlation lags / = where I depends on the statistical 
properties of the random processes («) and c rfprf (n) (namely for / > Ithe variance in the 
estimate even asymptotically does not attain the Cramer-Rao lower bound, a well known 
concept for the person skilled in the art, as it is shown hereby in said reference: B. Porat, 
"Digital processing of random signals", Englewood Cliffs, NJ, Prentice-Hall, 1994. 

In another preferred non-limitative embodiment, the following weighted sum of ML 
estimates can be considered: 

Pml =2><A>a( / ) (22) 

where K represents the maximum correlation lag at which the cross and autocorrelation (17) 
and (18) are estimated. Similar to the estimator (21) the optimal weighting coefficients 
y, can be computed based on the accuracy (expressed in terms of mean square error) of each 
ML estimate fi M (/) . In practice y, will increase with the accuracy of the estimate (/) , 
rendering the contribution of the most accurate estimates to the sum (21) more significant 
than that of less accurate estimates. Similar to estimator (21) the efficiency of the 
estimates fi M (l) will decrease as the lag / increases. Thus also in this estimator a maximum 
lag Zcan be determined for the estimate to be asymptotically efficient (see the said 

reference: B. Porat, "Digital processing of random signals ", Englewood Cliffs, NJ, 
Prentice-Hall, 1994) where we remark that it would be in general K*L. The derivation of 
the optimal weight coefficients a, and b, , or y, although of theoretical interest, results in far 
too complex a matter to allow a simple closed form expression suitable for a realistic 
implementation and is hence beyond the scope of this invention. However, it can be shown 
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that since the accuracy of the estimates p ML (l) as well as the accuracy of the estimates 
<f> dc (l) and ^(/) generally decreases as the correlation lag / increases the coefficients y 9 and 
a j and b t will decrease accordingly. 

In the light of the above considerations the least complex ML estimator with the 
maximum estimate accuracy is that of equation (20) with / = 1 , i.e.: 



n*=m 

Notice that estimator (23) has also the advantage of minimizing the delay needed to 
compute the cross and the autocorrelation at numerator and denominator respectively, 
reducing the complexity and the need of additional memory storage. 

1 0 4.3.2. Practical Means for Improving and Tracking the Beamforming Factor Estimate 
There might be practical circumstances where the beamforming factor p can stay 
constant over several slot periods. This may occur when the BS optimizes the beamforming 
weight vector w on the basis of the path angles 9 P and the average energies of the paths. In 
such cases since those quantities can stay approximately constant over several slot periods 

15 the BS will not modify the beamforming weight vector w so that the factor p p = w H a(^) for 
/>=!,..., J> will stay approximately constant during the same time period. This fact can be 
used to improve the accuracy of the estimates (19) — (23) by averaging/filtering the 
estimates of the cross- and auto-correlation at fixed lag / * 0 (17) and (18) computed over 
successive slot periods, namely over the period where the factor p is constant. Then for that 

20 averaging-filtering operation to be properly implemented the UE has to be informed 

whenever the beamforming factor^ is abruptly modified either because the angles of the 
paths have changed and/or the BS has modified the beamforming weight vector w , in order 
for the UE to reset/reinitialize the averaging-filters memory. Indeed if the averaging/filter 
memory is not reset, the UE would account for cross-correlation and autocorrelation 

25 estimates (16) — (17) corresponding to an out-of-date beamforming factor p . However, one 
may notice that there might be no need to reset the averaging-filter memory in the case of 
slow variations of the beamforming factor p . Indeed appropriate filtering of successive 
cross and auto-correlation estimates (17) and (18) inherently performs the tracking of slow 
variations of the beamforming factor £ . For this purpose we suggest the implementation of a 
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low complexity first-order low pass ICR filter with a suited forgetting factor (first-order IIR 
filtering is well known by the person of ordinary skill in the art). Eventually, in order for the 
BS to notify the UE about abrupt changes in the transmit beamforming, a signaling 
overhead is necessary though. If the standard under consideration, e.g. UMTS, does not 

5 envisage to support such signaling from the BS to the UE, either the averaging filters have 
to be removed or the UE needs to implement an additional mechanism able to detect 
whether an abrupt change occurred in order to make a decision on whether the averaging 
filter memory is to be reset. Such an additional detection mechanism as well as the cross and 
auto correlation averaging-filters, clearly come at the expense of complexity increase that 

10 might impact the UE design. However, if due to complexity limitation the averaging-filters 
are not implemented, the BS will not have to signal abrupt beamforming factor changes to 
the UE avoiding additional network overhead. The UE shall then update the estimate of the 
beamforming factor every slot, independently of the previous slot estimates. Furthermore 
one may observe that in this case the UE can keep estimating the factor J3 even 

15 independently of the presence of transmit beamforming. Indeed when no transmit 

beamforming is performed at the BS, ft will reduce to a real factor representing the transmit >. 
amplitude offset between the CPICH and DPCH. 

4.3.3. DPDCH Channel Estimation at Different DPCH Symbol Rates 

So far, for the sake of simplicity in the derivation we have assumed at a particular instance 

20 the spreading factor of the DPCH to be equal to 256, in order for the DPCH and the CPICH 
to have the same symbol rate (i.e. the same number of chips per symbol). However, as 
mentioned previously the spreading factor of the DPCH can take on values in the set {2, 4, 
8, 16, 32, 64, 128, 256, 512} so as the DPCH symbol rate generally differs from that of the 
CPICH whose spreading factor is always fixed equal to 256. Hence in a preferred non 

25 limitative-embodiment, the instantaneous ML channel estimates (8)-(10) generalize as 
follows: 

CcptchW = s %tch( n )ycptchW = c(n) + v^in) 

^(0 = ^(0^(0 = wVW)+v^(0 
where, given the DPCH spreading factor M , the index i corresponds to the i -th DPCH 
symbol spread over M chips, and n corresponds to the n -th CPICH symbol always spread 
30 over 256 chips. Then since to any CPICH symbol corresponds 256/ M DPCH symbols in 

order for the previous algorithm to be used either the CPICH or the DPCH ML channel have 
to be interpolated to match the symbol rate of the other channel estimates. 
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Moreover, in deriving the ML instantaneous channel estimates (8) — (10), both the 
CPICH and the DPCH channels are implicitly assumed constant over the longer symbol 
period of the two channels, namely the Doppler spread is always assumed to be significantly 
lower than the lowest symbol rate allowed by the standard concerned. We shall briefly show 
5 that this assumption holds in practice by considering the case of the UMTS standard. Indeed 
the lowest symbol rate arises when the spreading factor is maximal, i.e. for M = 512 . In the 
UMTS standard the chip rate is fixed equal to 3.84MHz (see: "3GPP Technical 
Specifications, Rel. 99, http://www.3gpp.or^specs/specs.html " incorporated herein by 
reference), so that the lowest symbol rate for the DPCH is equal to 3.84MHz/512= 7.5 kHz. 

10 The UMTS standard also envisages supporting wireless mobile services where the UE 
velocity can be up to 500km/h, i.e. in the presence of a Doppler spread at most equal to 
1kHz that is still largely below 7.5kHz for the above assumption to be valid. Then in a 
preferred non-limitative embodiment, interpolation is done to adapt the ML instantaneous 
CPICH and DPCH estimates to the lowest symbol rate between said CPICH and DPCH. In 

1 5 particular, when the DPCH spreading factor M is smaller than 256, the interpolation to the 
CPICH symbol rate is obtained by averaging DPCH ML static channel estimates over a 
number of symbols equal to 256/Af, whereas when M=512, being the DPCH symbol rate 
lower than that of the CPICH, the average is taken on the CPICH over two consecutive ML 
instantaneous channel estimates. One can notice that this embodiment has the advantage of 

20 not coloring the estimation noise, i.e. noise samples after interpolation are still uncorrelated, 
so that equation (16) still holds. 

4.4. DPDCH Channel Estimation Algorithms 

In this section we summarize the steps of the joint CPICH-DPCH MAP DPDCH channel 
estimation algorithm as well as the joint CPICH-DPCH low-complexity DPDCH channel 
25 estimation algorithm derived in sections 4.2 and 4.3 respectively. 

4.4.1. Optimal Joint DPCH-CPICH MAP DPDCH Channel Estimation Algorithm 

In this section we provide a summary of the main steps of implementing the optima 
joint DPCH-CPICH MAP DPDCH channel estimation algorithm described above. The 
following steps are to be performed to provide an optimal MAP DPDCH estimate according 
30 to (11). 

1 . Compute the ML instantaneous CPICH and DPCH channel estimates c epfch (n) and 
(0 according to equations (8)— (10) and the generalization (23). 
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2. Perform interpolation of the ML instantaneous DPCH and CPICH channel estimates 
obtained above to the lowest symbol rate of said DPCH and CPICH, i.e.: 

• if the DPCH spreading factor M < 256 the DPCH ML channel estimates are to be 
interpolated in order to match the CPICH symbol rate, 

5 • if the DPCH spreading factor Af = 512 the CPICH ML channel estimates are to be 

interpolated in order to match the DPCH symbol rate, 

• if the DPCH spreading factor M = 256 no interpolation is needed. 

The ML instantaneous channel estimates thus interpolated are collected in the 
following vectors: 

10 K P icH=lc ep icH(n-K epich +\) ... c cpicA (n)] T 

^dpch^\c dpch (i-K dpch ^r\) ... C dpch (i)f 

To further reduce the complexity we limit the estimates c dpch (n - K dpch + 1),..., c dpch (n) to 
the estimates over the DPCCH c dpccb (i-K dpcch +i),...,c, pccA (0 in order to avoid the 
implementation of a DD mechanism necessary to provide pilot symbol estimates 
1 5 over the DPDCH. The vector c m = [c T dpch c T epich ] T is formed. 

3 . Build the optimal linear prediction filter f according to equation (13) either by 
implementing one of the techniques described in section 4.3 

4. Filter the ML instantaneous channel estimates e m = [c^ c T cptdl f obtained at step 3 
with the linear prediction filter computed at step 4 according to equation (12) to 

20 obtain a MAP DPDCH channel estimate c dpdch _ MAP {k) 

5. If the DPCH spreading factor M < 256 , the above DPDCH channel estimate 
CjpjcH-MApik) is to be interpolated at the DPCH symbol rate to be properly used. The 
least complexity interpolation technique consists of repeating the estimate 
Cdpdd>-mp( k ) 256/M times, if M <> 256 and averaging over two successive estimates 

25 Cdpdch-MAP (*) if M - 5 12 . In a general preferred embodiment more sophisticated 

interpolation techniques can be adopted (e.g., linear, quadratic, spline interpolation) 

6. Repeat steps 1 to 6 for all paths p » l,...,p assigned to the RAKE receiver fingers. 

4.4.2. Joint CPICH-DPCH Low Complexity DPDCH Channel Estimation Algorithm 
In this section we summarize the main steps of the low complexity joint CPICH-DPCH ML 
30 optimal DPDCH channel estimation algorithm. The following steps are to be implemented 
to estimate the DPDCH channel c^Ofc) at time k within the UE. 
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1 . Compute the ML instantaneous CPICH and DPCH channel estimates c^ fch (ri) and 
c dpch (i) according to equations (8)-(10) and the generalization (23). 

2. Perform interpolation of the above obtained ML instantaneous DPCH and CPICH 
channel estimates obtained above to the lowest symbol rate of said DPCH and 
CPICH, i.e. 

• if the DPCH spreading factor M < 256 the DPCH ML channel estimates are to be 
interpolated in order to match the CPICH symbol rate, 

• if the DPCH spreading factor M = 512 the CPICH ML channel estimates are to be 
interpolated in order to match the DPCH symbol rate, 

• if the DPCH spreading factor M = 256 no interpolation is needed. 

3. The ML instantaneous channel estimates thus interpolated are collected in the 
following vectors: 

*cpick = tfcptch (w - K cptch - 1 + 1) c cp ich (*)f 

*dpch=lc dpch (n-K dpch +\) ... c dpeA (n)] T 
with / > 0 . In our implementation we shall set / = 1 and = K dpch . To further reduce 
the complexity we limit the estimates c dpch (n - K dpch + 1),..., («) to the estimates over 
the DPCCH c dpcch (n - K dpcch + 1),..., c dpcdt {n) in order to avoid the implementation of a DD 
mechanism necessary to provide pilot symbol estimates over the DPDCH. 

4. Build CPICH auto correlation and DPCH CPICH cross correlation at lag / > 0 of the 
form (17) and (18), i.e., 

and 

n»ro 

5. Build an ML estimate of lie beamforming complex factor >£as 

Z C ^W^("-0 

where we shall set 1 = 1 according to equation (23). 
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6. Build an accurate (e.g. MAP) CPICH channel estimate c epich {k) of the true CPICH 
channel c^ik) by suitable low-pass filtering the ML CPICH instantaneous channel 
estimates c cp/cA =[c c ^(«-^ teA -/+l) ... c cpich (n)f . To reduce the complexity of the 
estimation filter we may consider only a part of the vector of the ML CPICH 

5 instantaneous channel estimates for low pass filtering. 

7. Build a DPDCH channel estimate c rfpdWi (*)of the true DPDCH channel as 

8. If the DPCH spreading factor M <256, the above DPDCH channel estimate 
Cdpdch( k ) = fiAa^cpich( k ) is to be interpolated at the DPCH symbol rate to be properly 

1 0 used. As previously mentioned, the least complexity interpolation technique consists 

of repeating the estimate c dpdeh {k) = p^c^k) 256/M times, if M <> 256 and averaging 

over two successive estimates c dpdch (k) = P^c^k) if M = 512 . In a general preferred 

embodiment more sophisticated interpolation techniques can be adopted (e.g., linear, 
quadratic, spline interpolation). 

15 9. Repeat steps 1 to 8 for all paths p = assigned to the RAKE receiver fingers 



In a non-limitative preferred embodiment, the CPICH multipath coefficient's estimate at 
step 6 is done by a linear prediction as described, e.g., in the document: J. Bsltersee, 
G. Fock, P. Schultz-Rittich, and H. Meyr, '^Performance analysis of phasor estimation 
20 algorithms for FDD-UMTS RAKE receiver", IEEE 6 th Symp. on Spread Spectrum 
Techonologies and Applications, NJTT, NJ, September 2000 (incorporated herein by 
reference). 

4.4.3. Relation to DPCCH Based DPDCH Channel Estimation 

Some authors (see e.g.: 'TSG RAN WG 4 meeting #17, R4-010594 Ericsson, Dedicated 
25 Pilots, May 2001," "TSG RAN WG 4 meeting #17, R4-010528 Nokia, Proposal for 
user-specific beamformingfor UTRA FDD, May 2001" incorporated herein by reference, 
and references therein, available at 

http://ww.3gpp.org/^ have addressed the problem of 

DPDCH channel estimation in the presence of transmit beamforming. Mainly they conclude 
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that since only the DPCH is affected by the transmit beamforming only the information 
available from the DPCH channel is to be used for DPDCH channel estimation under those 
circumstances. In particular when DD DPDCH ML channel estimates cannot be used due to 
complexity limitations, they recommend performing DPDCH channel estimation only 
5 exploiting the a-priori knowledge of the DPCCH pilot symbols. We shall notice that when it 
is viable both the algorithms described here by making use of information available from 
both the CPICH and the DPCCH have two major advantages with respect to a DPDCH 
channel estimation only based on the DPCCH pilots. First, the scarceness of pilot symbols 
provided by the DPCCH (from 1 to a maximum of 16 depending on the slot format see the 

10 hereby incorporated reference "3GPP Technical Specifications, Rel. 99, 

http://www.3gpp.org/sp ecs/snecs.html" V limits the accuracy of DPDCH channel estimate 
because of lack of sufficient noise and interference suppression. Secondly, as mentioned 
above, within a DPCH time-slot the DPCCH pilot symbols are not transmitted continuously 
as the CPICH pilot symbols, but they are time-multiplexed with the unknown data symbols 

15 comprising the DPDCH. In the presence of high Doppler spread, i.e. at high UE velocity, 
such structure of the DPCH time-slot might pose the problem of estimating the channel path 
coefficient c^ ch (A) during the absence of known pilot symbols as during the DPDCH period, 
unless the CPICH pilot symbols which are continuously transmitted are not exploited as 
well. Indeed if the Doppler spread is very low compared with the slot rate, the channel 

20 c dpdch (*) will vary so slowly within a slot period that it could be considered approximately 
constant over that period. Thus the channel estimate provided over the DPCCH can be 
considered valid for the whole DPDCH period. In other words there is no practical need of 
updating the DPDCH channel estimate at a higher rate than the slot rate. Conversely, if the 
Doppler spread is of the same order of magnitude as the slot rate or even larger, the 

25 assumption that the DPDCH channel is approximately constant over a slot period is no 
longer valid. For instance, referring to the UMTS standard, a slot comprises 2560 chip 
periods yielding a slot rate of 1.5kHz whereas, under extreme cases, the system is supposed 
to cope with Doppler spread up to 1kHz. Therefore, in such cases the channel estimated 
over the DPCCH cannot be considered a valid approximation of the true channel c dpdch (k) for 

30 the entire successive DPDCH period. Under these circumstances one should resort to 

interpolation and prediction techniques or yet implement a DD mechanism s imil ar to the 
one previously described that uses hard decision taken on the DPDCH data symbols as 
known pilots to update the channel estimates in order to update at the appropriate rate the 
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channel estimate over the DPDCH period between two consecutive DPCCH bursts of 
known pilot symbols. A DPDCH channel estimator, including such a tracking mechanism 
for the channel estimates is likely to be too computationally demanding to be feasible. 
Furthermore, under certain practical circumstances, such as when the DPCCH logical 
5 channel consists of only a single pilot symbol per slot, that approach would definitely yield 
poor performance. The joint CPICH-DPCH MAP DPDCH channel estimation algorithm 
represents a solution to this problem since it implicitly performs prediction of the true 
DPDCH channel based on the whole available information adapting the prediction filter to 
the actual Doppler spread and number of DPCCH known pilot symbols. Similarly, also the 

1 0 second low complexity algorithm proposed here inherently exploits the presence of the 
known pilot symbols over both the DPCCH and CPICH logical channels and allows 
tracking of the DPDCH channel variations even at high Doppler spread. The latter algorithm 
results are more attractive from a practical implementation standpoint though due to the low 
complexity involved and the nearly optimal performance. Indeed the second DPDCH 

1 5 channel estimation algorithm, even in its least complexity instance, builds an estimate j} of ) 
the beamforming complex factor /?, jointly accounting for the DPCCH and CPICH available • 
information (i.e. all available known pilot symbols), and provides up-to-date estimates of i: 
the DPDCH channel by multiplying the estimated factor 0 by the channel estimate c cptcA (n) 
computed relying on the CPICH pilot symbols. The continuous provision of known pilot 

20 symbols over the CPICH channel can be easily exploited by the UE to track the variations 
due to the Doppler spread of the channel path coefficient c(k) by continuously updating the 
CPICH channel estimate (k) . The tracking of the DPDCH channel is inherently 
obtained as the product c dpdeh (k) = Aa**«(*) . Finally we shall remark that within our second 
algorithm the DPCCH pilots are not directly used to compute a DPDCH channel estimate 

25 but only to compute an estimate of the factor p. If necessary and/or if the UE can afford the 
increased complexity hard decision taken, the DPDCH can be used as pilot symbols to 
improve the estimate of P as represented by the equations derived above. 

4.4.4. Algorithms Working Limitations 

In order to setup a communication in any mobile communication system, a mobile terminal 
30 UE at power-up needs to listen to the network and acquire a certain number of parameters 
related to the operator, the BS with which the UE will start the communication, the type of 
service required etc. Prior TO the useful data transmission the UE has primarily to perform a 
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certain synchronization procedures to acquire a BS. Within the general framework of DS- 
CDMA systems suitable for implementing a RAKE as receiver at the mobile terminal UE, 
there exists a preliminary synchronization phase that provides along with other parameters 
and time reference estimates, also estimates of the time-delays r p , and the average energies 
5 associated with the paths comprising the channel impulse response from the BS concerned. 
In the UMTS standard the mobile terminal UE generally estimates the delays of both paths 
and average energies by exploiting the CPICH. Based on those estimates the RAKE receiver 
fingers can be properly setup, as it is well known by the person skilled in the art. According 
to the methods of the invention described above, we have implicitly assumed that when the 

10 receiver can detect a path and estimate the associated delay on the CPICH during the 

synchronization mechanism, a corresponding path at the same delay exists on the DPCH 
even though with a different complex amplitude. As we have seen, this assumption allowed 
the derivation of a DPDCH channel estimation algorithm with low complexity, good 
performances, and high flexibility, able to cope with the problem of DPDCH channel 

1 5 estimation in the presence of transmit beamforming. 

In the light of the above it will be noted that under certain circumstances the energy 
associated with a path might be high over the CPICH channel and fairly low over the DPCH * 
channel. In particular a path can be strong enough to be detected over the CPICH channel by 
the synchronization mechanism, while being too weak to be detected over the DPCH 

20 channel. Conversely, it may happen that a path that carries significant useful signal energy 
on the DPCH channel is actually undetectable on the CPICH channel, and therefore not 
detectable by the synchronization mechanism operating based only on the CPICH channel. 
The first case would lead to an estimate of J3 with a very small magnitude. Thus, the 
algorithms described above would still work properly and the RAKE receiver would simply 

25 drop the finger as the corresponding DPCH channel path's energy becomes too small. The 
second case is much more critical because there would be no way to exploit the information 
provided by synchronization mechanism operating based only on the CPICH channel. 
Under these circumstances any of the two DPDCH channel estimation algorithms would 
properly operate and one should first search for paths associated with the DPCH channel by 

30 implementing a synchronization mechanism based on the DPCCH. Such a further 

synchronization mechanism would involve a significant complexity increase along with 
several real-time management complications, and would likely lead to poor performance 
due to the limited number of DPCCH pilot symbols per slot as we have previously seen. 
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5. DPDCH Channel Estimator Architecture 

In this section we provide architecture schematics for the implementation of the previously 
described DPDCH channel estimation algorithms. 

5.1. Optimal Joint DPCH-CPICH MAP DPDCH Channel Estimation Algorithm 
5 Architecture 

Fig.3 represents a schematic architecture for an implementation of a RAKE receiver 
comprising the proposed DPDCH channel estimator using the algorithm summarized by 
steps 1-5 described in section 0. In particular Fig.3 describes a portion of a single finger of 
the RAKE receiver, comprising: 
10 • A delay recovering block denoted DELAYJRECOVER compensating the delay 
associated with the path assigned to the finger of interest 

• A descrambling block denoted DESCRAMB 

• A CPICH despreading block, associated with the CPICH channel denoted 
CPICHDESPREAD 

15 • A DPCH despreading block associated with the DPCH channel denoted 
DPCHJDESPREAD 

• A CPICH symbol derotation block, denoted CPICHJDEROT, removing the modulation 
by the n-th CPICH pilot symbol 

• A demultiplexer block denoted DEMUX to recover according to the current DPCH slot 
20 format, the pilot symbols comprising the DPCCH and separate the data symbols 

comprising the DPDCH to be decoded 

• A DPCCH symbol derotation block, denoted DPCCH JDEROT, removing the 
modulation by the z-th DPCCH pilot symbol 

• Two interpolation blocks, one for the CPICH denoted INTERP1 and another for the 
25 DPCH denoted INTERP2. The first interpolation block INTERP1 is enabled when the 

DPCH spreading factor M =512 . The second interpolation block INTERP2 is enabled 
only when the DPCH spreading factor M < 256 . No interpolation block is enabled when 
M = 256 (see step 2 section 4.4.1). 

• A joint CPICH-DPCH MAP DPDCH channel estimation low pass filter denoted 

30 DPDCH_MAP_CH_EST, involving a delay of L CPICH symbol periods, and providing 
the estimates c 4-toft . M4P (*) according to step 4 of section 4.4.1. 
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• A delay line of L(256/M) DPCH symbol periods denoted DELAY necessary to 
temporally align the MAP DPDCH channel estimate c dpdch _ MAP (k) with the DPDCH data 
symbols 

• An interpolation block denoted INTERP3 to adapt the estimate at the CPICH symbol 
5 rate to the DPCH symbol rate, according to step 5 of section 4.4. 1 . 

• A block denoted CMPLX_CNJ, computing the complex conjugate of the MAP DPDCH 
channel estimate after interpolation to the DPCH symbol rate at the output of the block 
INTERP3, necessary to the RAKE combiner, as it is well known by the person skilled in 
the art 

10 • A multiplier block denoted DPDCH_CHEST_DATA_MULT that computes the product 
of the time-aligned DPDCH symbols and the corresponding complex conjugate DPDCH 
channel estimate. The product result is then sent to the RAKE receiver combiner, 
denoted RAKE_COMBINER. 

15 The scheme depicted in Fig.3 is to be replicated for all RAKE fingers. The output of the 
multiplier block DPDCH_CHEST_D ATAJVTULT of each finger is summed up by the 
RAKE combiner RAKE_COMBINER to form the output of the RAKE receiver as is well 
known by the ordinary skilled person. 

In practice, in a preferred embodiment, nearest neighbor interpolation will be adopted 

20 for the implementation of all the interpolation blocks INTERP1, INTERP2, and INTERP3, 
as it is shown to be a proper solution for this purpose in terms of complexity and 
performance trade-off. Nearest neighbor interpolation is the least complexity interpolation 
method. It actually saves complexity at the expense of a higher signal distortion than that 
resulting with more sophisticated interpolation techniques. However it can be shown that the 

25 actual amount of signal distortion yielded by that interpolation method is significantly 

below the noise and interference for practical ranges of signal, interference, and noise levels. 
Therefore, the additional signal distortion due to nearest neighbor interpolation is actually 
negligible compared with the effects of the noise and interference so as not affecting the 
receiver performances. 

30 5.2. Joint CPICH-DPCH Low Complexity DPDCH Channel Estimation Algorithm 
Architecture 

Fig.4 represents a schematic architecture for an implementation of a RAKE receiver 
comprising the proposed DPDCH channel estimator using the algorithm summarized by 
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steps 1-8 described in section 4.4.2. In particular Fig.4 describes a portion of a singlefinger 
of the RAKE receiver, comprising: 

• A delay recovering block denoted DELAY_RECOVER compensating the delay 
associated with the path assigned to the finger of interest 

5 • A descrambling block denoted DESCRAMB 

• A CPICH despreading block, associated with the CPICH channel denoted 
CPICHJDESPREAD 

• A DPCH despreading block associated with the DPCH channel denoted 
DPCH_DESPREAD 

10 • A CPICH symbol denotation block, denoted CPICH JDEROT, removing the modulation 
by the w-th CPICH pilot symbol 

• A demultiplexer block denoted DEMUX to recover according to the current DPCH slot 
format, the pilot symbols comprising the DPCCH and separate the data symbols 
comprising the DPDCH to be decoded 

15 • A DPCCH symbol denotation block, denoted DPCCH JDEROT, removing the 
modulation by the z-th DPCCH pilot symbol 

• A CPICH channel estimator low pass filter denoted CPICH_CH_EST, involving a delay 
of L CPICH symbol periods, and providing the estimates c cpfch (it) according to step 6 of 
section 4.4.2. 

20 • The actual estimator of the beamforming factor J3 9 which comprises two interpolation 
blocks, one for the CPICH denoted BSfTERPl and another for the DPCH denoted 
INTERP2. The first interpolation block INTERP1 is enabled when the DPCH spreading 
factor M = 512. The second interpolation block INTERP2 is enabled only when the 
DPCH spreading factor M < 256 . No interpolation block is enabled when M = 256 see 

25 step 2 section 4.4.2. 

• A delay line of L CPICH symbol periods denoted DELAY1 necessary to temporally 
align the CPICH channel estimate c^ik) with the up-to-date estimate p of the 
beamforming factor fi 

• A delay line of L(256/M) DPCH symbol periods denoted DELAY2 necessary to 

30 temporally align the CPICH channel estimate (k) with the DPDCH data symbols 

• A low-pass-filter denoted LPF1 for the DPCH-CPICH cross-correlation estimate and an 
LPF denoted LPF2 for the CPICH autocorrelation estimate to be enabled when the UE is 
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aware that the beamforming factor fi stays constant. The LPF memories are to be reset 
as the UE detects that the beamfoiming complex factor has abruptly changed (either via 
signaling from the BS or via a suitable detection mechanism implemented within the 
UE). Both filters LPF1 and LPF2 are not to be used if the UE has no means to detect 
5 abrupt changes of the beamforming complex factor. 

• A multiplier denoted BETA_CPICH_CH_EST_MULT computing the product 

PML C cplch 

• An interpolation block denoted INTERP3 to adapt the estimate at the CPICH symbol 
rate c dpdeb (k) = p^c^lc) obtained at step 7 to the DPCH symbol rate, according to step 8 

10 of section 4.4.2. 

• A block denoted CMPLX_CNJ, computing the complex conjugate of the DPDCH 
channel estimate after interpolation to the DPCH symbol rate at the output of the block 
INTERP3, necessary for the RAKE combiner, as is well known by the person skilled in 
the art 

15 • A multiplier block denoted DPDCH_CHEST_DATA_MULT that computes the product 
of the time-aligned DPDCH symbols and the corresponding complex conjugate DPDCH 
channel estimate. The product result is then sent to the RAKE receiver combiner, 
denoted RAKE_COMBINER. 

20 The scheme depicted in Fig.4 is to be replicated for all RAKE fingers. The output of the 
multiplier block DPDCH_CHEST_DATA_MULT of each finger is summed up by the 
RAKE combiner RAKE_COMBINER to form the output of the RAKE receiver as it is well 
known by the ordinary skilled person. 

In practice, in a preferred embodiment, also in this case nearest neighbor interpolation 

25 will be adopted for the implementation of all the interpolation blocks INTERP1 , INTERP2, 
and INTERP3, as is shown to be a proper solution for this purpose in terms of complexity 
and performance trade-off. Indeed the same considerations on nearest neighbor 
interpolation had done in section 0 on the joint CPICH-DPCH MAP DPDCH channel 
estimation algorithm hold in the case of the joint CPICH-DPCH low complexity DPDCH 

30 channel estimation algorithm as well. 
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6. Concluding Remarks 

It is to be understood that the present invention is not limited to the aforementioned 
embodiments and variations and modifications may be made without departing from the 
spirit and scope of the invention as defined in the appended claims. In this respect, the 
5 following concluding remarks are made. 

It is to be understood that the present invention is not limited to the aforementioned 
UMTS application. It can be used within any application of DS-CDMA where the 
communication standard envisages the presence of two logical channels, i.e. a common pilot 
channel, comprising symbols or signal waveform(s) known at the mobile terminal UE, and a 

10 dedicated channel comprising unknown data symbols and, in a preferred but not limitative 
embodiment, also dedicated known pilot symbols or signal wavefonn(s). 

It has to be understood that the present invention can be applied in the uplink 
communication, from the mobile terminal UE to the base station BS, when the mobile 
terminal UE the employs transmit beamforming. 

15 It is to be understood that the methods according to the present invention are not 

limited to the aforementioned implementation. There are numerous ways of implementing 
functions of the methods according to the invention by means of items of hardware or' 
software, or both, provided that a single item of hardware or software can carry out several 
functions. It does not exclude an assembly of items of hardware or software or both carrying 

20 out a function. For instance, algorithm steps can be combined, thus forming a single 

function without modifying the method of channel estimate with transmits beamforming in 
accordance with the invention. 

Said hardware or software items can be implemented in several maimers, such as by 
means of wired electronic circuits or by means of an integrated circuit that is suitably 

25 programmed, respectively. The integrated circuit is to be part of a receiver that can be 
contained in a computer, in a mobile communication handset, in a base-station, or in any 
other communication system device. Said receiver comprises means of the receiver adapted 
to make all necessary operations to support a certain type of communication services (i.e. 
operations of de-scrambling, de-spreading, etc.) said means being hardware or software 

30 items as stated above. In a similar fashion, the first estimator within the receiver comprises 
maximum likelihood estimation means adapted to make the ML instantaneous estimates of 
steps 1 and 3, interpolation means adapted to make rate matching as in step 2, optimal MAP 
DPDCH channel multipath coefficient estimation means adapted to make the DPDCH 
estimates of step 4, interpolation means adapted to make rate matching as in step 5, as 
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previously described, said means being hardware or software items as stated above. In the 
same way, the second estimator within the receiver comprises maximum likelihood 
estimation means adapted to make the ML instantaneous estimates of steps 1 and 3, 
interpolation means adapted to make rate matching as in step 2, correlation estimation 
5 means adapted to make the correlations estimates of step 4, beamfonning complex factor 
estimation means that are adapted to make the factor complex estimate of step 5, a first 
multipath coefficient estimation means adapted to make said CPICH multipath coefficient 
estimate of step 6, a second multipath coefficient estimation means that are adapted to make 
said DPDCH multipath coefficient estimate of step 7, interpolation means adapted to make 

10 rate matching as in step 8, as described previously, said means being hardware or software 
items as stated above stated. Of course, as mentioned above, said means can be an assembly 
of items of hardware or software or both carrying out a function, or a single item carrying 
out several functions. 

The integrated circuit comprises a set of instructions. Thus, said set of instructions 

15 contained, for instance, in a computer programming memory or in a stand-alone receiver 
memory may cause the computer or the receiver to carry <5ut the different steps of the 
estimate methods. 

The set of instructions may be loaded into the programming memory by reading a 
data carrier such as, for instance, a disk. A service provider can also make the set of 
20 instructions available via a communication network such as, for example, the Internet. 

Any reference sign in the following claims should not be construed as limiting the 
claim. It will be obvious that the use of the verb "to comprise" and its conjugations does not 
exclude the presence of any other steps or elements besides those defined in any claim. The 
article "a" or "an" preceding an element or step does not exclude the presence of a plurality 
25 of such elements or steps. 



